Incipient magnetism in the cubic perovskites MgCNi 3 and YBRh 3 : A comparison 
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Using density-functional-based methods, we have studied the effects of incipient magnetism in 
the cubic perovskites MgCNiz and YBRhz. Our results show that (i) at the equilibrium volume, 
both MgCNiz and YBRhz alloys remain paramagnetic and (ii) at expanded volumes, only YBRhz 
shows the possibility of a ferromagnetic phase with a local magnetic moment larger than 0.25 /is 
per Rh atom. 
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The coexistence of magnetism and superconductivity, 
predicted to be incompatible |]J due to pair-breaking ef- 
fects of magnetic scattering of electrons Q], has been ob- 
served in a wide variety of materials such as the cuprates 
SOL the ruthenates |BE3|- the heavy fermion compounds 
[3 111, the organic superconductors 0, 0], the boro- 
carhide s fill Il2l ] and the transition-metal intermetallics 
ZrZn 2 [3f 

One of the recently discovered superconducting mate- 
rial, the cubic perovskite MgCNiz [Tj| with a supercon- 
ducting transition temperature Tc = 8K, shows subtle 
signatures of incipient magnetism in the form of spin- 
fluctuations. The density of states of MgCNiz has a 
strong van-Hove singularity, primarily composed of Ni 
3d states, just below the Fermi energy 
suggests that the material may be close to a magnetic 
instability. Using the rigid band approach, one expects 
~ 0.5 holes to drive MgCNiz to ferromagnetism [ltill7j . 
However, experiments find that hole-doped MgCNiz, 
with the hole-doping achieved via Fe or Co substitutions 
in the Ni sub-lattice d EEO, El E3| , and vacancies 
[23, E3 or B substitutions in the C sub-lattice [2^|, re- 
mains non-magnetic. 

A variety of transition-metal-rich cubic perovskite ma- 
terials have been synthesized [13] and tested for super- 
conductivity. Confirmed report of superconductivity ap- 
pears only for YB x Rhz, for which resistivity character- 
ization finds Tc of ~ IK jSjJ. Magnetic measurements 
show that YBRhj, remains a Pauli paramagnet over a 
wide temperature range, 4.2 < T < 300 K 0J. The den- 
sity of states of YBRhz reveals that the density of states 
at Ep, N(Ep), is dominated by the Rh 4d states, with 
little or no contribution of the B 2p states [23, |23| ■ For 
expanded volumes, the Ef is expected to lower in energy. 
Such a movement in Ep is expected to enhance N(Ep), 
thus raising the possibility of satisfying the Stoner cri- 
teria. At the equilibrium volume, model calculations for 
MgCNiz show that the Stoner criteria is far from be- 
ing satisfied [la 13 Cj , but is found large enough to induce 
spin-fluctuations. 

Given that MgCNiz and YBRhz are perovskite su- 
perconductors, and that incipient magnetism resides in 
MgCNiz in the form of spin-fluctuations USE!, it is in- 
teresting to compare the propensity of magnetism in the 
two alloys as a function of volume. The present study 
of ordered MgCNiz and YBRhz alloys is a step in that 
direction. 



The unpolarized, spin polarized, and fixed-spin mo- 
ment calculations for ordered MgCNiz and YBRhz al- 
loys are carried out, self-consistently, using the Korringa- 
Kohn-Rostoker method in the atomic sphere approxima- 
tion (KKR-ASA). The calculations are scalar- relativistic 
with the partial waves expanded up to l max — 3 in- 
side the atomic spheres. The exchange-correlation ef- 
fects are taken into account using the Perdew and Wang 
parametrization j^. The core states have been recal- 
culated after each iteration. The overlap volume result- 
ing from the blow up of the muffin-tin spheres was less 
than 15%. The integration of the Green's function over 
the energy to evaluate moments of the density of elec- 
tronic states was carried out along a semi circular con- 
tour comprising of 20 points in the complex plane. For 
the Brillouin zone integrations, the special k— point tech- 
nique was employed with 1771 k— points spread in the 
irreducible wedge of the cubic Brillouin zone. The con- 
vergence in the charge density was achieved so that the 
root-mean square of moments of the occupied partial den- 
sity of states becomes smaller than 10~ 6 . 

The calculated equilibrium lattice constants for 
MgCNiz and YBRhz are found to be 7.138 a.u. and 
7.919 a. u., respectively. For MgCNiz the calculated 
lattice constant is consistent with the previous first- 
principles reports [13], but the value is underestimated by 
1% when compared with the experiments [l3 |. Surpris- 
ingly, for YBRhz the calculated lattice constant appears 
to be overestimated in comparison with the experiments 
[27l I33I which is not so common. A previous the- 
oretical estimate of the lattice constant of YBRhz also 
finds an overestimated value l2Sl. Such a discrepancy 
also exists for ZnCNiz alloys j^H Si- For the latter, 
it is argued that the material ZnCNiz subjected to ex- 
periments may be off-stoichiometric in C sub-lattice |3fi| . 
Thus, it can be inferred that the material YBRhz, which 
is reported to be a superconductor j^J, may be deficient 
in B content, and if so, following the case of MgCNiz one 
may expect higher Tc for the stoichiometrically ordered 
materials. 

The densities of states of MgCNiz and YBRhz, calcu- 
lated at the respective equilibrium lattice constants, are 
shown in FigQ] For both the materials, the density of 
states at Ep is dominated by the transition-metal states, 
Ni in MgCNiz and Rh in YBRhz. For MgCNiz, the 
van-Hove singularity is well reproduced, and it is about 
0.06 mRy below Ep, in agreement with the earlier band- 
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Figure 1: The total and sub-lattice resolved partial densities 
of states of MgCNis (upper panel) and YBRI13 (lower panel) 
, calculated at their equilibrium lattice constants. The thick 
line represents the total density of states while dotted, dashed 
and thin lines represent partial densities of states for Ni (Rh), 
Mg (Y), and C (B). The inset is a blow up of the states near 
the Fermi energy which is shown as a vertical line through 
the energy zero. 



structure reports [ifj- For YBRh 3 , the van-Hove singu- 
larity falls at Ep. In the itinerant model of magnetism, 
a high N(Ep) is essential for magnetism, and thus it ap- 
pears that both MgCNis and YBRh^ alloys are close 
to ferromagnetic instability. In the case of YBRhz the 
total energies of both spin polarized and unpolarized 
calculations at the equilibrium lattice constant remain 
degenerate, making YBRh% into a paramagnetic mate- 
rial. The calculated paramagnetic N(Ep) for MgCNi^ 
and YBRI13 are 14.562 and 11.212 state/ Ry — atom re- 
spectively, which are consistent with the previous reports 
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The characteristic difference in the densities of states 
of MgCNia and YBRh$ is the absence of non-metal 2p 
states at Ep in YBRh^. This compares well with its iso- 
electronic counterpart MgCCoz- However, MgCCo 3 is 
found to be a weak ferromagnet with Co having a local 
magnetic moment of 0.33 \xb |3sL Ifliij. Next, we describe 
the effects of volume expansion on the incipient mag- 
netism in MgCNi 3 and YBRh 3 . 

For MgCNis, up to a 15% increase in lattice constant 
from the equilibrium showed no local magnetic moment 
at the Ni site. The results of our calculations are sum- 
marized in Figl2l where we show the density of states as 
a function of lattice constant in MgCNi 3 over an energy 
interval of -0.2 < E < 0.2. Though N(E F ) increases 
steadily with increasing lattice constant, the peak, which 
is characteristic of the Ni 3d bands, does not change its 
position on the energy scale, consistent with the observa- 
tions made by Rosner et. al. [IF| . However, the increase 
in N(E F ) is not sufficient to induce a magnetic phase 
transition in MgCNis even at expanded volumes. 

For YBRfi3, our calculations reveal a sharp ferromag- 
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Figure 2: The total density of states of MgCNi^ for various 
lattice constants (as indicated in the left inset), calculated as 
described in the text. The density of states is shown over an 
energy interval —0.2 < E < 0.2 with respect to the Fermi 
energy. The vertical line through the energy zero represents 
the Fermi energy. The total density of states at Ef as a 
function of lattice constant is shown in the right inset. 
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Figure 3: The total magnetic moment as well as its contri- 
butions from the Y, B and Rh sub-lattices in Y BRhz as a 
function of lattice constant, calculated as described in the 
text. 



netic transition with respect to expanded lattice con- 
stants. In FigEl we show the total magnetization as a 
function of lattice constant for YBRh^. The individual 
contributions to the magnetization from Y, B and Rh 
sub-lattices in YBRI13 are also shown in FigEI We find a 
sudden appearance of local magnetic moment larger than 
0.25 [is per Rh atom for lattice constants larger than 15% 
of the equilibrium value in YBRh^. Small magnetic mo- 
ments, which are anti-ferromagnetically coupled to the 
magnetic moment at the Rh site, also develop at the Y 
and B sites due to the hybridization with the 4d states 
of Rh. 

As a cross check to the self-consistent calculations, 
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Figure 4: The Magnetic energy AE(M) (in mRy) as a func- 
tion of magnetization M (in [ib) for MgCNi^ (upper panel) 
and YBRhz (lower panel), calculated as described in the text. 



as described above, we have carried out fixed-spin mo- 
ment calculations j4Q| to estimate the magnetic energy 
AE{M) =E{M) - E(0), where E(Q) and E(M) are the 
total energies for the paramagnetic phase and ferromag- 
netic phase with a magnetic moment M respectively, for 
both MgCNi 3 and YBRh 3 at various lattice constants. 



In Fig 31 we show the calculated magnetic energy 
AE(M) as a function of M for various lattice constants 
for both MgCNi 3 and YBRh 3 . From FigH it is clear 
that MgCNiz remains paramagnetic throughout, as no 
minimum, other than for M=0, appears in AE(M) vs. 
M curve. However, for YBRhz a second minimum ap- 
pears in the AE(M) vs. M curve, indicating a ferromag- 
netic transition. Hence, both self-consistent and fixed- 
spin moment methods show a greater propensity of mag- 
netism in YBRhz than in MgCNi^. The existence of 
an enhanced magnetic behavior in YBRh^is due to a 
weak hybridization of the metal d states with that of 
the non-metal p states in comparison with MgCNi^. In 
MgCNi3, one finds that the C 2p states and the Ni 3d 
states hybridize strongly with significant contribution to 
N(Ef) from the C sub-lattice. 



In conclusion, first-principles, density-functional-based 
calculations find YBRh^ to be closer to a magnetic in- 
stability than MgCNiz with respect to lattice constant. 
The enhanced propensity of magnetism in YBRh^ may 
be responsible for lowering the value of Tq, in spite 
of high density of states at Ep. The superconductor 
YBRhz thus falls in the same class as the refractory 
compounds VN, where spin-fluctuations lower the Tc 
drastically. 
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